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ABSTRACT 


The NPS computer simulation model was used to investi- 
gate the sputtering by 1.0 keV argon bombardment of clean 
and oxygen reacted monocrystalline titanium, vanadium and 
niobium. Variations in yield, energy and angular distri- 
butions and in particular, multimer formation, were studied 
as a function of oxygen density and location. Simulation 
results show a significant decline in substrate yield as 
the oxygen coverage is increased, regardless of location. 
Further, there is a marked preference for multimer forma- 
tion by lattice fragmentation, rather than by recombination, 
for all three metals. The percentages of multimers formed 
by fragmentation were found to increase with increasing oxy- 
gen density for all three metals. Previous results on face- 
centered cubic crystals concerning the yield per layer and 
the effects of channeling were confirmed for body-centered 
cubic and hexagonal closed-packed crystals. For compari- 
son, Simulation runs were repeated at ion energies of 0.5 
and 2.0 kev for a titanium with oxygen atoms located in 


three-fold sites with a C(2x2) density. 
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I. BACKGROUND 


A. OVERVIEW 

Sputtering is the ejection of atoms from a target by 
bombardment with an energetic projectile ion. The incoming 
ion collides with atoms of the solid, thereby losing a 
portion of its energy and momentum. If the energy gained 
by a struck atom exceeds the binding energy holding it in 
position in the lattice a primary recoil atom will be 
created. This recoil atom will then collide with other atoms 
in the material, distributing its energy via a collision 
cascade. An atom at or near the surface will be sputtered 
if the energy associated with the normal component of its 
velocity exceeds the surface binding energy of the material. 

The phenomenon of sputtering can be subdivided into two 
principle categories. Transmission sputtering is the ejec- 
tion of atoms from the rear of a thin target following pene- 
tration by the ion. Sufficient energy must be transported 
through the target to allow atoms to overcome their binding 
energy and escape. Back-sputtering is the ejection of atoms 
from the front of the target and depends upon sufficient 
energy being deposited in the surface layers to allow 
ejection. This thesis will be concerned with the latter. 

First discovered by Grove in 1853 and Faraday in 1854 
as the deposition of metal atoms on the glass walls of a 


gas discharge tube, it was not until 1902 that Goldstein 











(Ref. 1] presented evidence that the sputtering effect was 
caused by positive ions of the discharge impacting on the 
metal cathode. Since its principle effect was to contaminate 
experimental environments and to erode apparatus, the majority 
of early work centered around ways to control or eliminate 
what was considered to be an undesirable side effect. The 
necessity of increasing the life expectancies of high voltage 
vacuum tubes and the recognition of the possible applications 
of sputtering to thin film coatings gradually spurred inter- 
est in understanding the mechanics involved. 

Early theoretical work was hampered by a paucity of 
experimental data. Reproducibility of results was also a 
major problem due to an incomplete understanding of the 
factors involved. Since the sputtering yield, Y (defined as 
the number of atoms ejected per incident ion), is critically 
sensitive to ambient pressure and surface contamination, as 
well as the flux density and angle of incidence of the beam, 
progress was dependent on the development and refinement of 


experimental techniques. 


B. EARLY EXPERIMENTAL WORK 

The effect of pressure on yield was first demonstrated 
by Penning and Moubis [Ref. 2]. They found that an increase 
in pressure increased the collision frequency between the 
ejected atoms and surrounding gases, resulting in the back- 
scattering of the escaping atoms back onto the target 


surface. By keeping background pressure below Gime TOLS, 
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they mere able to obtain reproducible results for ion ener- 
gies in excess of 500 eV. Below 500 eV, the scouring action 
of the incident beam was insufficient to maintain a clean 
surface. 

The elimination of adsorbed gases and surface oxide 
layers as a critical factor in determining the true yield 
was described by Arifov, et al., in 1963 [Ref. 3]. Yonts 
and Harrison also presented evidence [Ref. 4] that surface 
recontamination from background gases was a significant 
factor in quantitative sputtering yield measurements. 
Criteria were therefore developed for conducting ‘clean' 
experiments. 

In addition to pressure and surface cleanliness, the 
following conditions must be met in order to obtain repro- 
ducible yields: 

(1) The ions must strike the target at a known angle 
of incidence. 
(2) The energy spread of the beam must be small. 
(3) All ions must be uniformly charged and mass separated. 
(4) On monocrystalline surfaces the lattice orientation 


must be adequately described. 


C. THEORETICAL DEVELOPMENT 

Pioneering work on the concept of an individual sputter- 
ing event on an atomic scale was conducted by Stark [Ref. 5]. 
Applying the conservation laws governing elastic collisions, 


he proposed two theories to explain known experimental 
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results. One, the "Hot-Spot" moaSt, considered the sputtered 
atoms to be the result of evaporation of target material 
from an atomically small region due to ion induced local 
heating. In his "Collision" model, sputtering events were 
the result of a series of binary collisions initiated by 

a single ion. His models were modestly successful in pre- 
dicting the energy distribution of the sputtered particles. 
Further refinements on the hot-spot model were attempted by 
von Hipple [Ref. 6], but only after the demonstration by 
Wehner [Ref. 7] of the effects of crystal structure on the 
yield did it become apparent that local heating alone could 
not explain the mechanisms involved in sputtering. 

Wehner's observation of spot patterns (angular distri- 
bution patterns characteristic of the crystallographic proper- 
ties of the substrate) obtained from sputtering of mono- 
crystalline targets revived interest in the collision approach. 
In 1952 Keywell [Ref. 8] made a first attempt to formulate 
Stark's collision model in terms of the neutron transport 
model developed for nuclear work. Harrison then developed 
a theory based on probability concepts as expressed by the 
idea of collision cross-sections [Ref. 9]. However, due to 
the large number of unknown parameters, these theories were 
limited in their quantitative application. 

As a further attempt to account for the angular distri- 
bution of ejecta, Silsbee [Ref. 10] proposed a focused 
collision model which allowed the transport of momentum in 


crystals along preferred directions. Available experimental 
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results showed that the yield from monocrystals depended 
sensitively on the crystallographic orientation of the inci- 
dent beam. This was explained in terms of the "holes" seen 
in crystal models when viewed from different planes and 

led to the development of various "Transparency" models. 
While focusing has been shown to be a contributing factor 

in cascade development within a crystal, it is inadequate 

to explain Wehner spots. 

As an alternative to the focusing collision theory, 
Lehmann and Sigmund [Ref. 10] proposed a model based on the 
Boltzman transport equations that required the target surface 
to have an ordered structure, but not long straight rows of 
atoms intersecting the surface as in the focusing model. 
Emphasis is placed on the role of the surface structure and 
in particular, surface binding energies in monocrystalline 
sputtering. As an ion strikes the surface a collision cas- 
cade is created via binary interactions. An atom in the 
cascade will sputter if its momentum component normal to 
the surface has sufficient energy associated with it to 
overcome the surface potential barrier. Thompson [Ref. 12] 
expanded on this idea and proposed that the surface attrac- 
tion for the escaping particle causes a refraction of its 
velocity vector away from the normal, resulting in a dis- 
torted angular distribution of ejected particles. Signifi- 
cant agreement with experimental results was found when 
projectile energies were much greater than the surface bind- 


ing energy. This approach also provided correct information 
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concerning both the angular and the energy distributions 
from polycrystalline targets. Thompson's [Ref. 12] work 
also incorporated the idea of focusons (momentum transport 
without mass displacement) into the theory. Robinson [Ref. 
13] has taken a similar approach while sates clin LOL 
channeling. 

While these theories have provided mathematical tools 
to aid experimentalists, they have been rather unsuccessful 
in predicting yields at low projectile energies and from 
Single crystal targets. Binary collision models have thus 
far been unable to account for the predominance of sputtered 
atoms from surface layers of the target. 

In general, these theories have attempted to reduce what 
is essentially a many-body, multiple interaction process 
into forms that are analytically tractable. An assortment 
of refinements followed, attempting to explain observed 
behavior. Full scale attack on the multiple interaction 
problem had to await the arrival of high speed computers with 
sufficient storage to handle the hundreds or thousands of 


atoms in a cascade. 


D. COMPUTER SIMULATION MODELS 

Paralleling the theoretical and experimental progress in 
the understanding of sputtering was the effort to model the 
basic mechanisms through the use of computer simulation. In 
1960 Gibson, Goland, Milgram and Vineyard [Ref. 14] developed 


a computer model to simulate the motion of atoms in a copper 
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monocrystal following the impact of an ion ona Single 
target atom. Atomic interactions were treated as binary 
collisions and the resulting cascades were developed by 
applying Newton's equations of motion. Ion penetration 
studies by Robinson and Oen [Ref. 15] using a similar 
program led to the discovery of channeling which was later 
experimentally verified. 

In 1967, based on earlier preliminary studies, Harrison, 
Levy, Johnson and Effron [Ref. 16] used a computer to simu- 
late the bombardment of a single copper crystal by argon. 

An important finding of this work was that the yield was 
almost entirely confined to the first three layers of the 

top surface. Subsequent refinements of the model [Refs. 17- 
20] ancluded the addition of potential functions that allowed 
for attractive, as well as repulsive interactions, creating 

a dynamically stable crystal. Harrison's results showed that 
focusons were not a major contributor to the yield. 

In 1978 Garrison, Winograd and Harrison [Ref. 21] pub- 
lished the results of a comprehensive study of atomic and 
molecular ejection from a copper crystal with adsorbed oxy- 
gen atoms. Their simulation of the system provided a de- 
tailed picture of the mechanisms of ejection and enabled 
them to determine whether molecules were ejected as a cluster 
from the surface or combined together in the space above the 
target. They were also able to determine the effects of 


adatom placement on molecule formation. Such results, when 
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coupled with laboratory data, should lead to a better under- 
standing of actual adatom locations on chemically reacted 
surfaces. 

While much work remains to be done, particularly with 
respect to constructing interatomic potentials that accurately 
represent the physical processes involved, computer simula- 
tion permits the study of sputtering phenomena with a preci- 
sion and flexibility that is impossible to duplicate in the 
laboratory. While it cannot replace actual experimentation, 
it can lead to a better understanding of experimental re- 


sults and provide a testing ground for theoretical concepts. 
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II. OBJECTIVES 


The primary objective of this thesis is to study the 
sputtering of titanium, vanadium and niobium monocrystals 
using the NPS computer simulation model, QRAD. In addition, 
changes in the yield and multimer formation are studied when 
the surface is reacted with oxygen. 

Yu [Refs. 22,23] has conducted comprehensive research 
on clean and oxygen reacted titanium, niobium and vanadium 
in an attempt to link several earlier theories on the 
mechanisms responsible for ionic multimer formation. In 
the Lattice Fragmentation Model first proposed by Benning- 
hoved [Ref. 24] molecular species (multimers) are emitted as 
complete fragments of the surface lattice during ion bombard- 
ment. In the Recombination Model [Refs. 18,25] only single 
atoms are ejected. Molecules are then formed through the 
attractive interaction between sputtered atoms from the 
same sputtering event, establishing their identities above 
the surface. The charge state of the molecules is not taken 
into account in these models. 

Yu has reported experimental evidence [Ref. 23] that the 
emission mechanism for a sputtered molecular species may be 
a function of its charge state. He proposed that for low 
oxygen coverages (< one monolayer) the formation of positively 


charged titanium and niobium molecules (MO, MO.) favored 


iL 


a — 
os 
<= 
i 
= @& 
= 
> 





the recombination model, while the formation of negatively 
charged species (MO , MO, and MO 3) followed the fragmenta- 
tion process in all three metals. In the case of vanadium 
it appeared that the positive oxide species were not formed 
exclusively by recombination, but by a combination of the 
two methods. 

Formation of positively charged ions was also found to 
be strongly dependent on oxygen coverage. When the number 
of oxygen atoms neighboring a metal atom increases, the 
frequency of sputtering an Mr and oh adjacent to each other 
naturally increases. Consequently the percentage of Mo* 
produced by fragmentation will increase. 

The formation of negative molecular ions by lattice 
fragmentation was found to be insensitive to oxygen coverage 
due to the kinetics involved in the formation proces. A 
molecule formed by recombination is usually left with excess 
energy. If this energy approaches the electron affinity 
of the molecule the negative ion will attempt to stabilize 
itself by ejecting an electron. If the excess energy is 
much greater than this, one approaches the dissociation 
energy, and the molecule will be prevented from forming at 
all. Thus the formation of a stable negative ion is most 
favorable when the components are close neighbors on the 
surface and eject with small relative velocities as in the 
lattice fragmentation process. Increasing the oxygen cover- 
age will only increase the probability of forming higher 


order oxide multimers. 
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Some discrepancies were noted in Yu's results for posi- 
tive ions and were attributed to the recombination of metal 
atoms and oxygen from adjacent sites. Thus the results were 
indistinguishable from fragmentation at higher oxygen 
coverages. All nine of the negative molecular ion species 
investigated by Yu showed strong evidence of formation by 
fragmentation. Since the energy balance for recombined 
negatively charged ions indicates that they tend to neutralize 
by electron ejection, this leaves fragmentation as the 
principal source of molecule formation. 

It should be noted that Yu's experimental procedure 
utilized an ion beam with a 45 degree angle of incidence to 
the polycrystalline surface. The mass spectrometer detector 
was also positioned at 45 degrees with respect to the surface, 
at right angles to the ion beam. Thus only a portion of the 
total yield was sampled. The data from the simulation studies 
that follow are based on the total yield from ion bombardment 
perpendicular to monocrystal surface. These data include 
neutral atoms, as well as ions, which are not collected by 
the mass spectrometer. 

While determination of the charge state of molecules is 
not possible with the computer model, the various mechanisms 
of ejection can be studied in detail to determine whether 
recombination or fragmentation is the major contributor to 
the multimer yield. Ejection times of the various constitu- 


ents are also available for correlation. 
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Additionally the energy and angular distributions of 
ejected particles, which can be related to experimental SIMS 
data, will be studied. Since the actual sites of adatom 
adsorption have yet to be determined, the effects of oxygen 
coverage and location on yield and multimer formation will 


also be investigated. 
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Ifift. MODEL DEVELOPMENT AND SIMULATION 


A. COMPUTER MODEL 

QRAD is a computer Steen program designed to track 
the motions of atoms in a microcrystal following ion bombard- 
ment. Trajectories Beep eotcanated by solving Hamilton's 
equations of motion using an average force method [Ref. 26]. 
It is a multiple interaction (MI) simulation following time 
step logic. Once the atomic masses, crystal characteristics, 
interatomic potentials and ion energy have been initialized, 
the program will develop and follow the resulting collision 
cascade until the energy of the atoms reaches a predetermined 
minimum. (See reference 24 for a comprehensive discussion 
of sputtering simulations.) 

Initially an ion of given mass iS projected toward a 
Surface impact point with a specified energy and angle of 
incidence. As the cascade develops, the positions and veloci- 
ties of the ion and lattice atoms are calculated at the end 
of each timestep. A timestep is determined by the velocity 
of the most energetic particle so that its displacement does 
not exceed 0.1 lattice unit ao) os To expedite calculations, 


particles are classified as either moving or non-moving de- 


pending on their energies and the forces exerted on them 


ithe lattice unit 1s a convenient method of expressing 
distances when dealing with different crystal types. An 
LU is equal to the lattice constant ‘a’ divided by root two 
for hep crystals and is equal to a/2 for bee crystals. 
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by neighbor atoms. Atoms which rise above the surface with 
sufficient momentum are reclassified and assigned to a 
tentative list of ejected atoms. A trajectory is terminated 
when the energy of the most energetic atom is such that the 
probability that either it, or other atoms, will escape the 
surface is Hae small. Final positions and velocities are 
then recorded. Multiple trajectories, simulating bombard- 
ment by an ion beam, can be run by computing trajectories 
for a set of impact points chosen to sample a representative 
area of the crystal surface as a function of its symmetry. 
Data on the sputtered atoms is then analyzed by separate 
programs (AN83 and PLOTAN) to determine multimer formation, 
ejection angle information, etc. 

As mentioned previously, multimers may result either 
from lattice fragmentation or by recombination above the 
surface. A multimer is considered formed when its constitu- 
ent atoms are ina favorable spatial position with suffi- 
ciently low relative kinetic energy to permit bonding. The 
total energy of the resulting system must be negative and 
is an indication of the group's binding stability. Absolute 
multimer yield is strongly dependent on the pair-potential 
in use and the range specified for atomic interactions 
(Ra), however the mechanisms for cluster formation have been 
found to be insensitive to the choice of the potential 
funeecion [Ref. 27]. 

Global results, suitable for comparison with experimental 


data, are obtained by averaging the yields from each run 
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over a number of trajectories. While the number of trajec- 
tories 1S quite small (100-200, depending on symmetry 
considerations), the results have been found to vary by 
less than 10% when compared with other independently chosen 
sets of impact points [Ref. 28]. These global results are 
used in determining the number of atoms ejected per single 
ion (ASI yield), energy and angular distributions, as well 
as the spatial distribution of yield from the impact Zone. 
The remainder of this chapter is devoted to the develop- 
ment of the potential functions and other parameters used 
by QRAD. Adatom nomenclature and placement are discussed 
in Section D. Trajectory sets and crystal impact points are 
covered in Sections F and G. The chemical properties of 
titanium, vanadium and niobium are covered in the following 


sections. 


B. TITANIUM 

Titanium is a dimorphic transition metal ranking ninth 
in abundance of the elements making up the earth's crust. 
Below about 880°C it exists as a Silvery gray hexagonal 
close-packed (hcp) crystalline material, while above 880°C 
it undergoes a phase transformation to a face-centered 
cubic (fcc) structure. The important physical and chemical 
properties of titanium are listed in Table l. 

Titanium reacts with oxygen to produce four well defined 
oxides, the monoxide TiO, the sesquioxide T1503, the dioxide 


Or titanic acid Ti0. and the trioxide or pertitanic acid 
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T1i0.. T150c, T1,0, and T1509, 5 have also been reported, but 
may in fact be metastable combinations of the others [Ref. 
Zo 

At least two distinct forms of oxygen coverage have been 
observed on T1(001) surfaces, 

1. Type 1 begins to form at low oxygen exposures of 
approximately 1.0 i and ls characterized by a low work 
function, very low electron stimulated desorption (ESD) 
and photon stimulated desorption (PSD) yields and a surface 
binding energy of 5.9 eV as determined by UV photoemission 
spectroscopy (UPS). Plateau, et al., had previously deter- 
mined this value to be 5.7 + 0.2 eV, also using UPS and 
ESCA [Ref. 38]. 

2. As the oxygen exposure is increased (~50 L) a second 
type of oxygen coverage is observed which has a higher work 
function, high PSD yield and additional binding energy in 
UPS. The effects of oxygen adsorption appear to saturate 
at about 70 L. 

Hanson [Ref. 37] concluded that type 1 may correspond 
to oxygen adatoms lying below the surface plane or in high 
coordination sites, possibly the three-fold hollows, which 
would account for the low desorption rate. 

Jonker [Ref. 39] has also observed two states of oxygen 
on the surface. The first is detected at low coverages and 


produces a well-ordered P(2x2) low energy diffraction 
2 -6 Pa -4 
1.0 L (Langmuir) = 10 Torr-sec = 1.33 x 10 Pa-sec. 
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(LEED) pattern. At higher oxygen exposures a second state 
emerges which has a disordered appearance. This was ex- 
plained by assuming that the sticking probability coefficient 
for oxygen is relatively constant up to a quarter monolayer 
coverage. At this point interference from previously 
adsorbed oxygens would cause it to change abruptly. The 
sticking coefficient would then remain constant for increas- 
ing exposures until the density of oxygen favored a change 
in the chemical bonding, resulting in the formation of a 
more stable metal oxide in the upper layers of the crystal. 
Kawasaki, et al. [Ref. 40], have also reached the conclu- 
sion that oxygen is first adsorbed on the surface at low 
doses, but as the exposure is increased, it is slowly ad- 
sorbed into the surface. This corresponds to the formation 
of an oxide, and is accompanied by changes in the crystal 
structure. Their results also indicate that at the lower 
exposures oxygen occupies either a doubly bonded A-top or 
a two-fold bridge position. The A-top position is also sup- 
ported by PSD experiments [Ref. 41] in which the oxygen 
atoms were observed to desorb normal to the surface, imply- 
ing that the Ti-O bonds are perpendicular to the surface. 
Experimental work by Singh, et al. [Ref. 43], confirms 
the formation of an oxide as saturation is reached, with 
oxygen atoms incorporated extensively into the metal lattice. 
(Here saturation is defined as the point where increasing 


the oxygen exposure no longer has an effect on the metal 
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crystal, i.e., the sticking probability is equal to zero.) 

Since the sputtering yield comes mainly from the top three 

layers, modeling of an oxygen saturated surface would have 

to be accomplished by using a titanium oxide crystal rather 
than titanium with oxygen adatoms. 

In view of the uncertainty of the actual locations of the 
oxygen atoms, a complete set of trajectories was run for the 
A-top, two-fold (bridge) and three-fold locations, in addi- 
tion to those on the clean surface. Figures 1 through 8 
present the various adatom locations on the metal crystal 


surface. (See Section D for a description of adatom notation.) 


C. VANADIUM AND NIOBIUM 

Vanadium and niobium are also transition metals, but 
unlike titanium, they have a body-centered cubic (bcc) 
structure. Their (110) planes, while not as close packed, 
do however have a hexagonal arrangement of atoms similar to 
the titanium (001) plane (Fig. 9). Their physical and 
chemical properties (Tables 2 and 3) are also quite similar 
to titanium, hence their sputtering mechanisms should bear 
some correlation to that of titanium. 

Haas, Jackson and Hooker [Ref. 44] in their LEED experi- 
ments on niobium, tantalum and vanadium have found that oxygen 
adsorption is the same on all three metals and is a function 
of oxygen exposure. At low oxygen coverages the oxygen 
patterns clearly indicate a six-fold symmetry characteris- 


tic of the substrate lattice. But, as in the case of 
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Fig. 5. Titanium with adatoms in P(2x2) three-fold "A" position. 
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Fig. 6. Titanium with adatoms in C(2x2) three-fold "A" position. 
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Fig. 7. Titanium with adatoms in P(2x2) three-fold "B" position. 
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Fig. 8. Titanium with adatoms in C(2x2) three-fold "B" position. Small 
dots indicate centers of atoms in the second layer. 


de) 














s in P(2x2) A-top position. 








COCOO 


&Y 


Y 
Q 


ye) 


Y 
Q 


yy 


g 
gy, 

& 

Y 


QOQDODG 


ADS 


COOCOC® 





TABLE 2 


Physical Properties of Vanadium and Parameters 


of Associated Potential Functions 


Atomic Number 23 

Atomic Weight 50.9415 

Melting Point 1s) SHC)p ae EO) 2 es 

Boiling Point Soguce e 

Valance +2 oa 

Crystal Type Dee 2 
2 3.0399. +,0.0003. 2 

Lattice Unit 1.51995 A 

BA AERIS 

Cohesive Energy ores a=) 


1°) 
Gas Phase Bond Lengths (A) 


V-V 260224 
Ti-O ives ee) 04 
O-O lie PAO pay 
Ar~O cea ie 


Gas Phase Bond Strengths (eV) 


V=V 2 oon 
v= 6.462 
O-O Deo4 5 
Ar~O 0.0049 


Surface Binding Energy 


V-O 5.52 eV 
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Gale. 


Calculated, 
see text. 





TABLE 3 


Physical Properties of Niobium and Parameters 
of Associated Potential Functions 


Atomic Number 41 REF 
Atomic Weight 92.9064 30 
Melting Point 2468 +10° C 30 
Boiling Point 4742° Cc 30 
Valance 2,3,4(2) ,5 30 
Crystal Type bee eter 
a 3.30656 + 0.00002 A 

Battic Unit ILE Sis Ae: A Caren 
R. 2.8636 A Calc. 
Cohesive Energy 7.5/7 eV 23 


7) 
Gas Phase Bond Lengths (A) 


Nb=-Nb 2.8584 30 
Nb-O 1.69 34 
0-0 ie 20752 Se 
Ar-O Speed 35 


Gas Phase Bond Strengths (eV) 


Nb-Nb Pe)! 34 
Nb-O Powe 32 
O=9 5.1545 35 
Ar-O 0.0049 30 


Surface Binding Energy 


Nb-O 6.67 eV Calculated, 
see text. 
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titanium, the reference index for the adatom mesh is uncer- 
tain. Two locations are possible on the vanadium and niobium 
(110) faces; the A-top and a rectangular four-fold site not 
present in titanium. In the four-fold position the adatoms 
could be placed either equidistant from the two pairs of 
neighbor atoms, or off to one side in a three-fold 'plus' 
position. Both would result in a high degree of bonding 
asymmetry. Since there is no corresponding site on the Ti 
(001) surface, the A-top position was chosen for comparison 
with titanium. 

No information is available on the binding energies of 
chemisorbed oxygen on vanadium or niobium. Since all three 
metals possess similar chemical properties, surface binding 
energies were calculated by scaling the gas phase M-O 
binding energies to the known SBE of oxygen on titanium. 
This resulted in estimated SBE's of 5.52 and 6.67 eV 


respectively. 


D. ADATOM COVERAGE 

Computations were run at two specific coverages, P(2x2) 
and C(2x2), for each of the adatom locations. These cover- 
ages correspond roughly to a quarter and a half monolayer 
of oxygen respectively. 

The designations for surface coverages were standardized 
im an article by E.A. Wood in 1964 [Ref. 45]. In this nota- 
tion the three dimensional unit cell is reduced to a diperiodic 


structure referred to as a net in which the unit area is the 
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unit mesh. The fourteen Bravais space lattices are thus 


reduced to five basic nets pictured below. 





Hexagonal 
Y = 120°, a=b 








a a a 
Oblique Rectangle Centered Rectangle 
heer 90 ,.a F b Y = 90°, a 7b Y = 90°, a # b 


The letter "P" is used to denote the primative mesh with 
equivalent points at the four corners, but none within the 
mesh itself. The letter "C" is used to identify a unit mesh 
that contains a central point. Numerals are used to indi- 
cate the repetition interval in the X and Z directions, 
referenced to the substrate unit mesh. Meshes rotated from 
the substrate principal directions are designated as above 
followed by the angle of rotation. Nac1(100) P(5x1)-45 


would represent a primative mesh located on the (100) 
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surface of NaCl, five units by one unit in area and oriented 
in the <ll0> direction. Meshes formed on other than sub- 
Strate atoms are further identified by their origin point, 
two-fold bridge, etc. (For simplicity the four number 
designation normally used for hcp crystals has been abbre- 
viated to three digits; eee (0001) = (001). The numbers 
moomoelves refer to the principal crystal directions (hkil) 
where 1= h + k.) 

In the case of the titanium (001) surface and the (110) 
surfaces of vanadium and niobium, the unit mesh is a 120 
degree rhombus. The unit mesh for the A-top positions is 
centered above the substrate atoms and is displaced for the 
two and three-fold positions. Figures l-ll depict the various 
coverages used, with the adatoms scaled according to their 
covalent radii. The two-fold locations are abbreviated "BR" 
(bridge) and the three-fold sites by "TF". Since there are 
two different three-fold sites on the hcp (001) surface, the 
additional identifier "TPA" or "TFB" is used to denote 
whether the adatom is fixed "above" an atom in the second 


layer of the substrate or centered over a channel. 


E. POTENTIAL FUNCTIONS 
1. General 
A variety of potential function forms are available 
in this program to describe the interatomic forces. The 
merits of each have been hotly debated in the literature, 


but Harrison has found [Ref. 42] that the global results are 
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not strongly dependent on the particular form in use. Indi- 
vidual parameters may be adjusted, if necessary, so that 
-absolute values match experimental results, but the basic 
mechanisms involved remain unchanged. The four potential 


functions used are as follows: 
a. Born-Mayer, which is strictly repulsive having the 


form 
Vir) = A exp(-br) ; 


b. Moliere, which is an approximation to the Thomas-Fermi 


screening function, and is also repulsive 


a) = (2, Z,e°/z) (0.35 exp (-0.3r/a) + 0.55 exp (-1.2r/a) 


+ 0.1 exp (-6.0r/a)} , 


where: 


LZ 


1/2) 2/3 
‘ 2 


a = 0.8853a,/ (2 + Z (Firsov radius) 


and 
2 II 
= h'/4n° me = (OR SS Ie (Bohr radius) 


The parameter "a" may be varied to produce a "modified" 


Moliere function. 
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c. Morse, which is both repulsive and attractive having 


the form 
V(r) = D, exp {-2a(r-r,) } - 2D, expt-a(r-r.)} , 


where re is the equilibrium separation distance between 
atoms, D, is the well depth and alpha is used to scale the 
distance in the potential function; and 
d. Composite Moliere-Morse, which consists of the two 

potential functions joined together smoothly by a cubic 
spline. 

The adjustable parameters can be related to the ion 
"size" and a collision "hardness", which in the case of the 
Born-Mayer potential correspond to the A and b constants. 
Thus one can tailor the frequency and hardness of collisions 
by appropriate selection of the parameters. At large separa- 
tions the Moliere and Born-Mayer functions have similar shapes. 
In all of the potentials the repulsive wall dominates the 
dynamics of collisions, while the attractive well determines 
whether an atom will escape the surface and be sputtered. 
Reference [42] contains a detailed discussion on potential 
function selection. 

2. Solid Phase Potentials 

The Ti-Ti interactions within the crystal were 
modeled using a standard Moliere-Morse spline (a = 0.10529 A) 
Lomaibicweter the anticipated hard collisions in the experi- 


o 
ment. R, was set at the lattice parameter value of 2.95 A 
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ad ee was adjusted to give the correct crystal cohesive 
eergy of 4.855 eV. The Morse alpha was set at 1.33 a 
s that the slope closely matched that of the Moliere. The 
tro were then connected by a spline between Ra = 0.68 and 
= om72 LU. Ras the distance at which the potential is 
truncated, was set at 2.20 LU to allow only interactions be- 
ween nearest and next-nearest neighbor atoms. Since the 
orces produced at greater separations are very weak, this 
election of R, allows calculation time to be minimized with 
egligible effect on the results. 

The Nb-Nb and V-V potentials were similarly chosen 


oi 


Co °o. 
nth alpha equal to 1.34A . and 1.62 A respectively. Ra 


fas set at 2.2 LU in both cases. R, was calculated as 43/2 
zimes. the lattice constant for the bcc structure. 

The metal (M)-Ar potential functions used for the 
oombardment of the crystal by argon were chosen as standard 
Moliere potentials since no binding is to be expected. Ra 
Was set at 1.7 LU for all runs. 

Morse potentials were used for the M-O surface 
interactions. Lacking any specific data in the literature, 
adatoms in the A-top position were placed at the M-O gas 


phase equilibrium separation distance above the substrate, 


although this is known to be too high. The well depth and 


alpha were then adjusted so that the correct surface binding 


energy was obtained. For the titanium two-fold bridge, gas 


phase equilibrium separation was Maintained between the oxygen 
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and its two neighbor titaniums which necessitated placing the 
oxygens 0.322 LU above the surface plane. The well depth 

was again adjusted for the correct SBE. Adatoms in the 
three-fold sites were placed in the center of the titanium 
triangles on the surface plane, which required increasing 

os ES a. 70 A. Ro for all was set just outside the nearest 
neighbor distance to prevent O-O interference. Parameter 
values for the various adatom positions are summarized in 
Table 4. 

The Ar-O potential used is a Born-Mayer with A = 14.67 
keV and b = ~-4.593 Ate This choice gave good results in 
previous work by Harrison, et al. [Ref. 21]. 

The O-O potential is another Morse with Dees Sila 
eV, ies imeos A and alpha = 2.691 al. values which were 


obtained from gas phase data and were also used by Harrison, 


et al. 
TABLE 4 
Solid Phase Potential Function Parameters 

ce] o.] ce] 

D, (ev) re (A) Alpha (A ~) Ra (A) 

Ti ATOP 5.90 We) 2 2.0 2.086 

Ti BRIDGE 29 5 1.62 20 2.086 
Die. Te le Sun 1. 70 2.0 2.086 

V ATOP Se 52 1.56 220 De yevs 

Nb ATOP 6.67 1.69 220 22810 
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3. Gas Phase Potentials 
The potentials used for interactions above the 
surface, as well as in determining multimer formation, are 
all of the pure Morse type. Parameters were selected from 
References 27, 32 and 33 and are summarized in Table 5. 
Appendix A (Figs. 15-22), contain plots of potential energy 


versus separation for most of the potential functions. 


TABLE 5 


Gas Phase Potential Function Parameters 


D, (ev) 2 (A) Alpha (ie) 
Ti-Ti 1.466 2.896 P7330 
Ti-O 5 tei (10 JEG) 410 2.000 
Nb-Nb 2G OU 25850 1.340 
Nb-O i605 1.690 2.000 
V-V ee oUt PL (3 OP 1.620 
V-O 6.462 13560 2.000 
O-O Siri) ses. eZ 2269 
MICROCRYSTAL 


Microcrystals composed of 676 atoms were used for each 
computer run. The atoms are placed by a lattice generator 
into the proper hep or bcc arrangement, resulting in a 
crystal 15 planes wide by 15 planes long by six planes deep 
(15x6x15). The first atom in the crystal is number two, 
one being the ion, and is located at (0,0,0,0). Atoms 2-114 


comprise the top layer, 115-226 layer two, 227-339 layer three, 
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340-451 layer four, 452-564 layer five and 565-676 layer 
six. Atom numbers above 676 correspond to oxygen adatoms. 
Crystal size selection is based on the problem of confine- 
ment which is the attempt to encompass the entire surface 
sensitive portion of the collision cascade within the target 
lattice. As the target size is increased the yield will 
asymptotically approach a maximum value. Past this point 
further increases in size will only affect computer run times. 
The 15x6xl5 crystal size was selected after making a series 
of runs on clean titanium using a high yield impact point. 
Yield versus ion energy and crystal size is depicted in 
Figures 12 and 13. As can be seen containment is effectively 


achieved for ion energies up to 3.0 keV. 


G. IMPACT AREA 

The impact area for each simulation is a rectangular 
array of 104 uniformly distributed impact points adapted from 
previous work [Ref. 46]. Figure 14 is a graphical represen- 
tation of the points and their coordinates. The actual length 
and height of the grid is automatically scaled by the pro- 
gram so that its dimensions match those of the interatomic 
Spacing in the X and Z directions. Thus a grid positioned 
on planes RBX = 6.0 and RBZ = 6.0 would have its four corners 
located on the centers of the atoms at (6,6), (8,6), (6,8) 
ana (o-c) « 

The impact points are slightly offset in the X and 2 
directions from true center to prevent repetitive sampling 


of symmetrically equivalent points. 
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Standard impact zone with 104 impact points. 
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maeeSaendent data sets can be generated by simply moving 
the origin of the grid. Selection of grid position is based 
on the surface symmetry resulting from adatom density and 
placement. Generally two sets of impact points are re- 
quired to sample a completely representative area of the 
target surface. Surfaces which could be mapped into one 
grid were still run with two sets, shifted by one plane, 


so that all runs consisted of 208 impact points. 


H. ION ENERGY 

All simulations were conducted with an ion energy equal 
to 1.0 keV. This energy produced a reasonable sputtering 
yield and allowed the runs to be conducted on a moderate 
sized crystal. Additional runs were conducted on the 
titanium C(2x2)TFB surface at 0.5 and 2.0 keV, for compari- 
son. This coverage produced the highest yield at 1.0 kev 
and therefore the differences in results as the ion energy 


1s varied should be more pronounced. 


I. PROGRAM MODIFICATIONS 
1. QRAD 
During the early stages of this investigation 
Harrison completed extensive modifications to QRAD resulting 
in improved energy conservation, timestep logic and output 
format. Considerable effort was expended proofing and de- 
bugging the new version. Additional code was also inserted 


to identify registered: pairs, metal-axygen pairs on the 
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undamaged surface, for later use in tracking multimer 
SL-igains. 
2. AN83 
This is the analyzer program used to determine vari- 
ous yield distributions and multimer data. Changes were 
made to make the program compatible with the new output 
format of QRAD. 
3. PLOTAN 
This program used the output from QRAD to generate 
graphic displays of the results. The program was modified 
to permit separate analyses of each atom type. It was also 
expanded to accommodate additional crystal types. The logic 
was modified to allow determination of the angular distributions 
from data generated by the non-symmetric sets of impact 
points in use. Plot scaling routines were also changed to 
handle the additional crystal types. 
4. ADPLOT 
ADPLOT is used to generate the pictorial representa- 
tion of adatoms and steps on the crystal surface. Modifica- 
tions were made to permit real time changes to the scaling 


of the plot size and to ddjust adatom radii. 
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IV. RESULTS 


A. GENERAL 

Detailed information on each run has been assembled into 
appendices following the main body of this thesis. Included 
are graphical displays of the ejection time distributions, 
atoms ejected per single ion (ASI), energy distributions, 
Spot patterns, and the ejection time differences between 
multimer components. Where applicable these have been broken 
down into separate substrate and adatom plots. Plot headings 
are generally self-explanatory, indicating ion energy, metal 
type and crystal face used for the particular run. For 
example “1.0 keV T1i(001)/Ar + Ox P(2x2)" indicates that the 
data pertains to bombardment of the 001 face of titanium 
by 1.0 keV argon ions. Ox P(2x2) denotes the presence of 
oxygen adatoms in a P(2x2) configuration. (M) or (OX) are 
used to identify whether the plot is for the substrate or 
for the adatoms in the particular configuration. The plots 
contain information based on a standard run of 208 impact 
points. On occasion data from a multiple set of 520 impact 
points is included to amplify significant characteristics. 
These are identified by the notation (IPX5). Combined plots 
of adatoms in the TFA and TFB positions are identified by 
TF (A+B). Discussions of the more noteworthy results of the 


experiments are contained in the remainder of this chapter. 
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Bee LleGD AND ASI 
The yields from the runs conducted on the clean crystal 
surfaces are summarized in Tables 6 and 7. All runs consisted 
of the standard set of 208 impact points with the exception 
of the expanded run on titanium which was composed of 520 
impact points. This larger run was used (ee test the sta- 
tistical validity of findings from the smaller sets. Of 
particular interest is the fact that the yield from the 
representative sample of 208 impact points for titanium 
is within 1.6% of the yield from a much larger run of 520 
impact points, demonstrating the validity of this rather small 
set. The significantly higher yield of argon from the bom- 
bardment of niobium and vanadium is to be expected in view 
of the lower mass ratio, 0.429 (versus 0.834 and 0.784 for 
Ti and V respectively) and the increase in the collision 
cross-section which is a function of the atomic number. 
The increased yield of V and Nb, as compared to Ti, may 

be explained by two characteristics: 

1. The more open nature of the bcc (110) face. While 
the overall cohesive energy is greater, the surface layer 
has fewer nearest neighbors and therefore fewer bonds must 
be broken for an atom to escape the surface. 

2. The fact that there are no channels present in the 
bec (110) face for bombardment at normal angles of incidence. 
Hence all the ion energy is deposited in the top layers of 


the crystal from which the vast majority of sputtered atoms 
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TABLE 6 


Yield and ASI (by Coverage) 


COVERAGE N % 
M (3) No (%) Nar Yu Yo 
CLEAN 
Ta (520) 945 > 10 1.8 - 
i (208) 374 = 2 1.8 - 
V 946 - 8 BS = 
Nb 682 = 52 SoS = 
P (2x2) 
Ti ATOP 291 (79) oe ay 7) 2 1.4 0.4 
V ATOP 616 (83) ZO) 6 3.0 0.6 
Nb ATOP 522 (82) 112 (18) 37 Dao O25 
Ti BR S520) 82 (20) 4 1.6 0.4 
ri TeA 324 (90) 38 (10) 2 6 OZ 
Ti TFB 386 (37) 58 (13) 4 1.9 Os 
C (2x2) 
Ti ATOP 2156062) 1-3 3-038) 2 L210 0.6 
V ATOP 457 (70) 198 (30) 6 ZeZ 0.9 
Nb ATOP 574 (70) 159 (36) 41 Es On7 
Ti Br 2/1 (65) 148 (35) 4 5.3 O45 7% 
Ti TFA 290 (77) S6. (23) 4 1.4 0.4 
Ti TFB(O.5) 233 (69) LOS (SL) 7 Pel 0.5 
Ti TFB(1.0) 368 (75) 124 (25) 6 1.8 0.6 
Pierre (2.0) 419 (78) 116 (22) 2 pare | 0.6 


NOTE: All ion energies are 1.0 keV unless otherwise noted, 
i.e., (0.5). The designations (208) and (520) 
refer to a single set of impact points and a 
combined group of 520 impact points respectively. 


Nar is the number of reflected argon ions. Nu 


and No are the actual counts of ejected atoms. Yu 


and Yo are the computed average yields, atoms/ion. 
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TABLE 7 


Yield and ASI (by Substrate) 


COVERAGE N % N % 
M (3) 0 (3) NaR Yu Yo 
CLEAN 
Tz (520) 945 - 10 eee ~ 
ji (208) By4 - 2 1.8 - 
mx 2X2) 
Tae ATOP Done 79 ) i 2 led 0.4 
Ti BR 331 (80) 82 (20) 4 1.6 0.4 
Ti TF <A+B> 355 (88) 48 (12) 3 1.8 Oe 
C22) 
Ti ATOP 215 (62) 133 “(38) 2. 10 0.6 
Ti BR 271 (65) 148 (35) 4 Les Oe? 
Ti TF <A+B> 329 (76) 105 (24) 5 1.6 0.5 
V CLEAN 946 = 8 4.5 - 
V P(2X2) ATOP 616 (83) 26 Cl) 6 3.0 OE=G6 
wer 2<2)) ATOP 457 (70) 198 (30) 6 2? 0.9 
Nb CLEAN 682 - 52 3.3 - 
Nb P(2X2) ATOP 522 (82) 112 (18) 37 D5 0.5 
Nisme (2X2) ATOP 374 (70) 159 (30) 4l 1.8 07 


NOTE: The quantity <A+B> for the three-fold sites 
refers to the combined average of the A and B 
positions of the adatoms. 


ye 





originate. The yield per layer data for adatoms in the ~A- 
top position is given in Table 8. Yield ratios for the 


other adatoms locations were found to be essentially identical. 


TABLE 8 


Sputtering Yield by Layer 


COVERAGE LAYER 1 LAYER 2 LAYER 3 LAYER 4 
CLEAN 
Ea (001) Sy) 0 0 0 
V (110) 900 44 8 0 
NB (110) 694 26 8 0 
P (2x2) 
ma (001) 285 5 iL 0 
ame cil O ) 566 41 7 2 
Nb (110) 490 24 8 0 
eC2x2 ) 
i (0 01) 213 2 0 0 
7 ee!) 429 24 8 1 
Nb (110) 347 24 3 0 


With respect to the various adatom coverages the metal 
yields for the A-top positions are fairly consistent at 
about 80% for the P(2x2) and 70% for the C(2x2). As the 
oxygen coverage is increased there are fewer direct impacts 
on the substrate atoms which begin to look more like a second 
layer to the incoming ions with a corresponding decrease in 
yield. This decrease is Summarized below where Delta Yu 
represents the drop in substrate yield as the coverage is 


increased from a quarter to a half monolayer. 
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TABLE 9 


Change in Metal Yield as a Function of Adatom Coverage 


COVERAGE DELTA Y¥,, (8) 
Ti P + CATOP -17 
V P + CATOP -13 
Nb P > C ATOP -12 
Ti P + CBR -15 
Ti P +> C TFA -13 
Ti P > C TFB -12 


As the adatom position is shifted from A-top to bridge 
to three-fold, the total yield (M+0) fluctuates with the 
maximum coming from the bridge positions. The metal yields 
are about the same percentage for both the A-top and bridge, 
but increase in the three-fold positions. The corresponding 
decrease in the fraction of oxygens ejected is reasonable 
in light of the added stability of the triply bonded positions. 

The increased metal yield in the bridge, and particularly 
in the TFB positions, can again be explained in terms of 
channel blockage. The bridge locations partially obstruct 
the channel in the hcp (001) plane, while the TFB position 
above the channel closes it to incoming ions. 

The Atoms per Single Ion graphs for selected cases are 
contained in Appendix B. The distributions for the substrate 
atoms are in general broader than for the adatoms. As the 
oxygen density is increased, there is little change in the 


substrate ASI distribution, but the maximum in the oxygen 


6 0 





distributions shifts toward an ASI of two. This is to be 
expected due to the larger number of oxygen atoms on the 
surface. The ASI distributions of the substrate broaden 
somewhat as the adatom positions are changed from the bridge 
to the three-fold, consistent with the increasing metal 


yield. 


°C. EJECTED PARTICLE ENERGY VERSUS COVERAGE 

It is interesting to note that as the adatom coverage 
in the Metop position is increased, the average energy of 
the ejected substrate atoms also increases. As noted above, 
this corresponds to a decrease in metal atom yield as well. 
Lower energy particles are apparently retained in the crystal, 
thus raising the average ejection energy. The most probable 
energy of the ejected metal atoms also is significantly 
higher when there are adatoms in the A-top position and in- 
creases with oxygen density as well (Table 10). This trend 
does not extend to the bridge and three-fold positions in 
titanium, but the most probable energies of the substrate 
atoms are more in line with those ejected from clean metal 
surfaces. 

There is a much closer correlation between the metal and 
oxygen average energies in the bridge and three-fold posi- 
tions. The multiple bonding involved appears to equalize 
the energies of the ejected particles to some degree. 

As the ion energy is increased there is a corresponding 


rise in the average ejection energy of the metal atoms, as 
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TABLE 10 


Average and Most Probable Ejection Energies 


COVERAGE 


CLEAN 


igel 
V 
Nb 


P (2x2) 


Ti 
ny 
Nb 


(ATOP) 


C (2x2) 


Ti 
V 
Nb 


(ATOP) 


TITANIUM 


P(2x2) BR 
C(2x2) BR 


P(2x2) TFA 
C(2x2) TFA 


P(2x2) TFB 
C(2x2) TFB 


<E 


Me (eV) 
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[Ey] MP (eV) 


<E,? (eV) 


10.062 
9.461 


On 997 
10.190 


8.055 
9.384 








TABLE 10 


Average and Most Probable Ejection Energies 


COVERAGE <Ey,? (eV) (Ey! MP (eV) <EQ? (eV) 
CLEAN 

jue 10.3 335 - 

V 8.9 4.0 - 

Nb 12 eS 4.5 - 
E22) (ATOP) 

aoe OF fan) ie. 

V CO 27 Tee: 11.4 

Nb INS aes} 9.0 eS. 
C(2x2:) (ATOP) 

Aube Ih ses Figs) 24 

V dee 8.0 TO 

Nb 15.9 9.0 38 
TITANIUM 
P(2x2) BR 103057 4.0 OrO Gz 
C(2x2) BR 9.815 S20 9.461 
P(2x2) TFA 9.852 S20 10.997 
@4Zx%2) TPA Ore t9:3 4.5 10.190 
P(2x2) TFB 107.236 3.0 CaO D5 

3.0 9.384 


CVzx2) TFB 974 3 
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TABLE ll 


Average Ejection Energy Versus Ion Energy 


ION ENERGY (keV) Ti (eV) Ox (eV) 
O5 8.059 225504 
a0 9.743 9.348 
20 13.146 20.316 


one would expect. The changes in the oxygen energy match 
the variation in yield (105, 124, 116). While the ratio of 
energies between Ti and O is essentially the same for 0.5 
and 2.0 keV, the oxygen energy is about double that for the 
1.0 keV ions at all coverages except the A-top C(2x2) where 
there is also an unexpected rise. 

The energy distribution curves are contained in Appendix 
C. The curves for the ejected oxygen atoms are typically 
flat for the bridge and the three-fold positions, showing 
no dominant ejection energy. However there is a pronounced 
Rayleigh-type shape to the curves for oxygen in the A-top 
position and they bear more of a resemblance to those of the 
substrate. This agrees with what was found earlier where 
the multiple bonds tend to cause a more even distribution of 


ejected particle energy. 


D. EJECTION TIME DISTRIBUTIONS 

Appendix D contains bar charts of the ejection time dis- 
tributions for representative runs. The distributions are 
similar for titanium and vanadium, but the average ejecticn 


time for niobium increases about 30 fsec to about 110 fsec 
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(1 femptosecond = 1071° segue, The distribution is also 


broader. There is essentially no change in any of the dis- 
tributions as oxygen is added to the surface, or increased 
from P(2x2) to C(2x2). Thus the ejection time appears to 
be mainly a function of the substrate mass and independent 
of adatom coverage and lattice type. 

The oxygen ejection time distributions show the same 
characteristics as the substrate atoms. The only significant 
difference was that oxygen ejection started approximately 
10 fsec later for titanium than for vanadium or niobium. 
Again the distribution was broader for niobium, reflecting 


the influence of the substrate. 


E. YIELD VERSUS ENERGY 

Since the C(2x2) TFB coverage produced the highest titanium 
sputtering yield, it was chosen to investigate the effects 
of changing ion energy. This coverage was therefore run 
at energies of 0.5 and 2.0 keV in addition to the 1.0 kev 
trial. The yield results and layer data are contained in 


Taores 12-and 13. 


TABLE 12 


Titanium Atoms Ejected and Yield 
as a Function of Ion Energy 


- Nu Yu No Yo Nar 
O25 233 dee 105 O25 
Ea 8) 368 Rees: 124 0.6 
2.0 419 220 116 0.6 
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TABLE 13 


Titanium Atom Yield per Layer 
as a Function of Ion Energy 


ENERGY LAYER L LAYER 2 LAYER 3 LAYER 4 
(kev) 
G25 232 I 
363 
20 398 16 1 


As can be seen the higher the energy, the higher the yield, 
but the relationship is definitely non-linear. With increasing 
energy the average ion penetrates further into the crystal. 
As a consequence a smaller percentage of its energy is deposited 
in the top layers. More energy is carried deep into the 
crystal where it is dissipated through collision cascades 
producing only lattice displacements and vibrations. 

The decrease in oxygen yield for the 2.0 kev ions is 
also probably caused by the deeper deposition of energy. 
The energy distribution plot for the ejected oxygens (Figs. 
54-56) clearly shows a decrease in the energy as compared to 
the 1.0 keV case. Hence fewer adatoms are receiving enough 


energy to surmount the surface binding well and escape. 


F. MULTIMER YIELD 
Tables 14 and 15 present the various multimer yields in 
terms of total numbers and as percentages. Tables 16 and 1/7 


contain a breakdown by species arranged by substrate type 
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TABLE 14 


Multimer Yield as a Function of Coverage 


TOTAL 
COVERAGE 1-MER 2-MER 3-MER 4-MER 5-MER MULTIMER 
Clean 
ma (520) 860 43 2 0 0 45. 
Ti (208) S50 10 2 0 0 eZ 
V 731 74 22 2 0 98 
Nb 585 48 13 2 0 63 
P (2x2) 
Ti ATOP 304 2S 4 0 0 29 
V <ATOP 529 77 2 5 0 94 
Nb ATOP 441 60 18 2 0 80 
Ti BR 350 26 3 0 0 29 
iy TRA Ba/ 10 i 0 0 all 
Ti TFB 389 ye} ME 2 0 26 
C(2x2) 
Ti ATOP 250 49 2 0 0 Syl 
V ATOP 416 99 3 3) 2 107 
Nb ATOP Si) 59 7 2 0 68 
Ti BR a35 yl 3 2 0 36 
T TRA 336 16 2 0 0 18 
T TB (.5) 310 14 0 0 0 jal 
Gy TPB (1) 437 18 5 0 0 23 
hapless. ( 2.) 490 Ly & 0 0 20 


NOTE: All ion energies are 1.0 keV except as indicated, 
ees. (0.5) . 
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TABLE 15 


Multimer Yield (Percentage) 


% YIELD 
THAT FORMED 
COVERAGE 1-MER 2-MER 3-MER 4-MER 5-MER MULTIMERS 
Clean 
ma (520) 90.4 9.0 0.6 0.0 0.0 9.6 
Ti (208) 97.9 348 1.6 0.0 0.0 6.9 
V 76.9 15.4 6.9 0.8 0.0 23 
Nb 80.5 bly raed See: oe 0.0 19.5 
P (2x2) 
Ti ATOP 83.2 LSS Se5 0.0 0.0 16.8 
V ATOP 14. 6 20.8 4.9 Peas | 0.0 28.4 
Nb ATOP 7.3 18.9 8.5 1.3 0.0 2S) 
Tc. BR Bioee2 12.6 Ze 0.0 0.0 aes 
Ta TRA 23 5.5 0.8 0.0 0.0 6.3 
Ti TFB 87.1 LOn4 (oe, 1.8 0.0 12.9 
ClZx2) 
Ti ATOP WO... 28.2 1.7 0.0 0.0 29.9 
V <ATOP Gp. 1 S0i.2 ibe! 1.8 i 34.9 
Nb ATOP 7230 22) 3.9 ie 0.0 245 
ra. BR B22 aes Data | 1.9 0.0 Seas 
A tye a by rs 89 .9 8.5 Ls 0.0 0.0 OE 
Ti TFB (.5) on. 7 8.3 0.0 0.0 0.0 8.3 
T? TRB (1.) 89.6 WS coe t 0.0 0.0 10.4 
Tie Tes (2. ) 93.6 6.4 jee 0.0 0.0 3h IL 


NOTE: All ion energies are 1.0 keV except as indicated, 
teers ee O« D) 
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TABLE 16 


Multimer Species Formed (by Substrate) 


COVERAGE M 
2 M, M, MO M0 M,0 MO. M0. M,0 
Clean 
ie (520) 43 2 0 - - - - = = 
Ts (208) 10 2 0 — = = = = a. 
P (2x2) 
Ti ATOP S 0 0 22 4 0 0 0 0 
Ti BR 8 0 0 18 3 0 0 0 0 
Ti TF <A+B> 7 0 0 8 at ip 0 0 0 
@iZ2x2) 
Ti ATOP 0 0 0 49 Z 0 0 0 0 
Ti BR 5 0 0 26 Z IL l 1 0 
Ti TF <A+B> 6 af 0 11 2 0 IL 0 0 
V CLEAN 74 ZZ 2 = - ~ ~ - - 
V P(2x2) 18 0 2 59 eZ 3 0 0 0 
te ( 252 ) 9 0 0 90 3 ie 0 Z 2 
Nb CLEAN 48 13 Z ~ - - ~ - = 
Nb P(2x2) g:: 6 0 42 TE, le 0 l 0 
NB C(2x2) 2 2 0 57 5 0 I 2 0 


NOTE: The quantity <A+B> refers to the average of the two 
three-fold positions. 
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TABLE 17 


Multimer Species Formed (by Coverage) 


RAGE 
COVERAG M. M. M, MO M,O M,0 MO. MoO, M,0. 
Clean 
Ti (520) 43 2 0 - - = = = = 
Ti (208) 10 2 0 = = = c= = _ 
V 74 22 2 - = = = = = 
Nb 48 183 2 - = = = ES = 
P (2x2) 
Ti ATOP 3 0 0 22 4 0 0 0 0 
V ATOP 18 0 2 59 i? 3} 0 0 0 
Nb ATOP 18 6 0 42 a i 0 1 0 
Ti BR 8 0 0 18 3 0 0 0 0 
Ti BFA 6 0 0 4 I 0 0 0 0 
Ti TFB Jah 0 0 12 af 2 0 0 0 
@i02x2) 
Ti ATOP 0 0 0 49 2 0 0 0 0 
V ATOP 9 0 0 90 & uf 0 2 2 
Nb ATOP 2 2 0 57 5 0 a 2 0 
Ti BR 5 0 0 26 2 iL iL rT 0 
Ti TFA 6 0 0 10 2 0 0 0 0 
toe ree (2.5) 4 0 0 10 0 0 0 0 0 
die Be (1. ) 7 IL 0 ll 2 0 2 0 0 
tr tre (2.) 9 0 0 8 2 0 l 0 0 


NOTE: All ion energies are 1.0 keV except as indicated, 
ieee, (O05) . 
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and adatom coverage. The multimer yield is seen to increase 
with oxygen coverage and is due primarily to the ejection 

of metal atoms with their attached adatoms. There is a good 
correspondence in the multimer yield between vanadium and 
niobium in the clean and P(2x2) coverage which then degener- 
ates as the coverage is increased to C(2x2). This can be 
attributed to the decreasing influence of the crystal lattice 
and increasing effect of the oxygen on the ejection mechan- 
isms. This will be covered in greater detail in the 


following section. 


G. MULTIMER ORIGINS 

One of the primary objectives was to study the origins 
of multimer components to compare fragmentation and recom- 
bination as modes of molecule formation. Table 18 presents 
data as to whether the constituent atoms were nearest-neighors 
in the crystal prior to ejection, or in the case of trimers 
or greater, whether a portion of the mdlecule was joined. 
In the table the letter 'C' represents the number of multimers 
whose components were coupled together on the surface; p-- 
partially connected multimers, N--those multimers whose 
eempenemtsewere not joined before ejection, tl=-the percentage 
that was formed by fragmentation and %2--the percentage 
that were formed either through total or partial fragmentation, 
where a portion of the multimer was joined on the surface. 


Table 19 contains similar data for trimers and higher order 
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TABLE 18 


Fragmentation versus Recombination 


(All multimers) 


COVERAGE c P N $i $2 
Clean 
Ti (520) 29 1 LS 64.4 66.6 
Ti (208) 6 1 5 S06 58.3 
V ak 7 20 2a 79.6 
Nb 4l L Pd 6521 66.7 
P (2x2) 
Ti ATOP 24 0 5 82.8 82.8 
V ATOP 76 6 12 80.9 87.2 
Nb ATOP 52 7 PAA 65.0 73.8 
Ti BR 23 2 4 79.3 86.2 
nL TEA 9 1 1 Sues 90.9 
aS EB oak 0 5 80.8 80.8 
G( 252) 
Ti ATOP 49 0 2 96.0 96.0 
V ATOP 86 2 19 80 .4 S252 
Nb ATOP 49 4 15 72.1 77.9 
Ti BR 29 2 5 80 .6 86.1 
Ti TFA ll 0 7 61.6 61.6 
Ti TRB (.5) 10 0 4 Jie 71.4 
Piemes {1.) 17 ae 5 T3329 7023 
Ti TRB (2.) is ie 6 65.0 10.0 


NOTE: C represents the number of multimers that were nearest 
neighbors on the surface prior to ejection; P, those 
partially connected and N, the number formed by 
recombination of non-adjacent atoms. $%1 is the 
percent formed by fragmentation, %2 the percent formed 
by partial fragmentation. All ion energies are 1.0 


keV except as indicated, i.e., (0.5). 


gil 








TABLE 19 


Fragmentation versus Recombination 
(Trimers and Greater) 


COVERAGE C P N $1 $2 
Clean 
Ti (52.0) 1 - 0 50.0 00 0 
Ti (208) 1 1 0 50.0 100.0 
Vi ny F, 0 70.8 100.0 
Nb 14 ils 0 93.3 MOO 0 
P (2x2) 
Ti ATOP 23 0 3 B38 5 88.5 
V ATOP sig 6 0 64.7 100.0 
Nb ATOP Le q) 0 65.0 100.0 
Ti BR 1 2 0 BS5 100.0 
Ti TFA 0 If 0 0.0 100.0 
Ti TFB S 0 0 100.0 100.0 
C22) 
Ti ATOP ap 0 IE 50.0 50.0 
V ATOP 6 2 0 T3528 100.0 
Nb ATOP 5 4 0 55.6 100 «0 
Ti BR 4 2 0 60.0 O00 
Ti TRA IL 0 il 50.0 50.0 
Tiere tl.) 3 1 ip 60.0 80.0 
is (2.) | 1 l B35 66.7 


NOTE: C represents the number of multimers that were 
nearest neighbors on the surface prior to ejec- 
tion; P, those partially connected and N, the 
number formed by recombination of non-adjacent 
atoms. %l is the percent formed by fragmenta- 
tion, %*2 the percent formed by partial fragmenta- 
tion. All ion energies are 1.0 keV except as 
micicated, 1.€., (0.5). 
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multimers. Table 20 presents the percentage of multimer 
species formed by fragmentation. 

The general trend is toward more fragmentation as the 
oxygen coverage is increased. This is due primarily to the 
fact that substrate atoms are ejected after being hit from 
underneath and carry off their attached adatoms as they leave. 
Many partial fragmentation molecules were formed when attached 
atoms exited the surface obliquely, knocking off oxygens from 
nieghboring atoms on the way out. This also appears to be 
one of the principal mechanisms involved in recombination 
of non-ajoining atoms. 

Figures 74-87 of Appendix E and Table 21 give the loca- 
tions and frequency of occurrence for all of the observed 
mechanisms involved in multimer formation. (In the table 
double entries refer to the a and b sections of the particu- 
lar drawing.) By far the greatest source of multimers is 
the dimer composed of a substrate atom and its attached 
adatom which is a fragmentation mechanism. Higher order 
multimers are a fairly rare occurrence, but over 65% are 
formed by fragmentation and virtually all involve at least 
partial fragmentation (Table 20). 

As the oxygen coverage is increased, the monoxide is 
still the prevalent species formed, but there is a shift in 
the mechanisms that make up the rest of the yield and new 
mechanisms are added. In the A-top and bridge positions 


the percentage of multimers resulting from fragmentation 
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TABLE 20 


Percentage of Multimer Species 
Formed by Fragmentation 


COVERAGE Mo M M MO MOS ee MO M.O 


3 4 Z 3 2 Za Se 
Clean 
ANGE 67 50 - - - - - - - 
V 74 7 oOo = - 2 = = = 
Nb 56 92 100 ~ - - - = - 
P (2x2) 
Ti ATOP 38 - = 82 100 _ = a = 
V ATOP 100 - 100 80 50 100 - - - 
Nb ATOP 83 100 = 57 50 yoo = 1P ” 
Ti BR 63 ~ - 89 Es) - - = _ 
fo A 100 - = 75 1P - - - - 
Ti TFB 64 = - 92 OOM LOo - ~ = 
Ti TF (A+B) 82 = - 84 50 100 - = - 
C1 2x2) 
Ti ATOP = - - 98 50 - - - - 
V ATOP 56 ~ - 81 67 1P - 100 100 
Nb ATOP 100 100 = 74 40 - X x - 
Ti BR 40 = - 92 100 100 1p 1p - 
Ta PRA a3 - ~ 80 50 = - = = 
Ti TFB 86 X = 73 100 = 50 - - 
Ti TR(A+B) 60 X - va 75 = 50 = ~ 


NOTE: A tack (-) indicates that the species was not formed, 
while "X" indicates that none were formed by frag- 
mentation. 1P indicates that only one of the particu- 
lar species was formed and that was through partial 
fragmentation. 
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TABLE 21 


Frequency of Occurrence of the Various 
Ejection Mechanisms 


COVERAGE iL 2 3 4 5 6 7 8 9 10 
Clean 
Ti (520) 2075) 91 2 973 1 1 - - - - 
Ti (208) BZ 2 2 1 = = - — 
V 87, 13 4 9 5 7 - a. 
Nb 19 8 10 4 6 2 8 - 2 - 
P(2x2) 
Ti ATOP cl ls 49 72 ~ - - - - 
V ATOP 6/0 6 1 = - 2 47/10 - - 4 
NB ATOP SAL. 43 3 3 4 - 24/7 Tel <5) 
Ti BR 4 AL l 1 1 Ae 2 deen 
Ti TFA AL 2 1 1k 5 1 = - - - 
Ti TFB 1 5 1 1 uf al 1 se ak = pill 
C(2x2) 
Ti ATOP J dhs 1 i = = - - - - 
V ATOP 5 l 9 al 2 = 74 Se es & 
Nb ATOP 2 - = - = 2 42 10 2 - 
Ti BR ik al 1 24 2 all ie askew Jit 
Ti TRA es 1 1 1” 6/2 1 l1- - 
tet) (.5) - - - 1 57 - yee 
fetes (1.) lee el Ae Fal0he 2 1 4/4 1/1 1a Racal 
Ti TFB (2.) l - 3/2 - - 5/1 = 1 - 


NOTE: The column headings refer to surface locations shown in 
Appendix E. Columns with double entries refer to the a and 
b sections of the individual figures. Ion energies are 1.0 
kev except as indicated, i.e., (.5). 
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increases, while it decreases for the three-fold positions. 
Those oxygens, sharing three bonds, tend to be left behind 
unless at least two of their nearest neighbors are also 
ejected. The greater oxygen density also presents a larger 
target to obliquely ejected fragments, increasing the 


probability of collision and subsequent capture. 


H. MULTIMER EJECTION TIMES 

Appendix F contains bar charts depicting the differences 
in ejection times between multimer components, both for 
metal-metal and metal-oxygen pairs. Essentially none of 
the pairs in any of the runs were ejected simultaneously, 
but all but about seven pairs were ejected within 100 fsec 
of each other and most within 50 fsec. This difference is 
about the same order of magnitude as the vibrational period 
of an equivalent diatomic molecule [Ref. 47], iio) == sec 
(100 fsec), and is not necessarily an indication that the 
atoms were ever disjoined at some point in time. 

Comparison of the charts for the different oxygen cover- 
ages reveals a correspondence between ejection time differ- 
ence and the percentage of multimers that were formed by 
fragmentation. As one would expect there is an increasing 


time difference between component ejection as the number of 


molecules formed by recombination increases. 


L. SPOT PATTERNS 
The spot patterns produced by ion bombardment of the 


crystal surfaces are contained in Appendix G. The ejection 


ies 








patterns correspond to the directions of the close-packed 
rows in the crystal, which for hcp are oriented at 307490 

and 150 degrees from the <001> direction. For bcc the direc- 
tions are 0 and 90 degrees. These lead to a six-fold symmetry 
pattern in hcp and a four-fold rectangular pattern in bcc. 
The patterns for vanadium and niobium are in quite good 
agreement with the expected results, but there is conewaee= 
able deviation in the case of titanium atoms. The spikes at 
90, 150 and 330 degrees are present, but there are extra 
Spikes at 0, 60, 180 and 240 degrees, while those at 30 and 
210 degrees are missing. A variety of different runs were 
conducted with additional impact point origins and permutated 
impact point sets, but all resulted in the same pattern. At 
present there is no apparent reason for this discrepancy. 
Second layer contributions would tend to broaden the spikes 
somewhat, but not change the overall appearance. 

The patterns for the metals with adatoms in the A-top 
position are essentially the same as for those from clean 
surfaces, even though they have considerably more energy as 
was pointed out earlier. There is some blurring, presumably 
due to adatom interference, but it is very minor. This 
blurring becomes more pronounced as the oxygen density is 
increased from P to Cc. The features of the metal spot 
patterns begin to fade when adatoms are placed in the bridge 
positions and are completely absent for oxygen in the three- 


fold sites, leaving a uniform angular distribution. 
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There is no marked tendency for oxygen to ejicee Per pela lcus 
lar to the surface in any of the locations, as was found by 
Hanson [Ref. 37]. The ejection angles were observed to in- 
crease slightly with increasing oxygen density, probably 
caused by the interference of other adatoms and increasing 
yield. This was also found to be the case for varying ion 
energy, where the average ejection angle reached a maximum 
for 1.0 keV bombardment corresponding to the peak in 


oxygen yield. 


J. EJECTION TIME VERSUS ENERGY 

The plots of ejection time versus energy are contained in 
Appendix H. As can been seen, the highest energy atoms are 
ejected very early in the process. As time progresses, the 
ejection energy falls rapidly with a few low energy stragglers 
coming off as the collision cascade dies and recrystalliza- 
tion begins. The distributions are independent of either 
adatom density or location. The duration of the ejection 
period is again seen to be a function of the crystal atomic 
mass as was found earlier in comparing ejection time dis- 
tributions. Increasing ion energy increases the energy of 
those atoms which are ejected early, but has no appreciable 
effect on the bulk of the ejected atoms. 

The energy distribution of ejected oxygens mimics that 
of the substrate, with the most energetic atoms coming off 
early in the cascade. This close correlation is most likely 


due to fragmentation where oxygens are ejected with their 
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neighboring metal atom(s). There is, though, a faster drop- 
off in energy than in the case of the substrate. Some 
broadening is observed with more energetic atoms being 
ejected later as the adsorption sites are shifted to those 


involving a higher degree of bonding. 


K. SURFACE LAYER EJECTION FREQUENCY 

Appendix I contains diagrams depicting the relative 
frequency of ejection of surface layer atoms. For simplicity, 
the plots contain information for a single set of 104 impact 
points. Also indicated are the locations of the impact 
zones that produced the ejections. 

The most obvious feature is the tendency for the damage 
area to extend vertically in the Z direction for the hcp 
crystals, while the pattern for bcc is much more symmetric. 
This can be explained in terms of the locations of close- 
packed rows. The two close-packed rows in bee are in the X 
and Z directions, however hcp only has one of three in the 
X direction, with the other two primarily in the Z direction. 

The plots also graphically demonstrate the decrease in 
substrate yield as the oxygen density is increased, as was 


mentioned in Section B. 


L. COMPARISON WITH EXPERIMENTAL DATA 
Tables 22 and 23 provide a comparison of experimental 
data on the ratios of multimer species formed by recombina- 


tion and fragmentation with that obtained by computer 
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TABLE 22 


Comparison of Experimental and Simulation Results 
on Multimer Formation Mechanisms for Titanium 


1% & P & P Cc P c 
RATIO YU (2x2) (2x2) (2x2) (2x2) (2x2) (2x2) (2x2) (2x2) 
ATOP ATOP BR BR TFAtB TFA+B ALL ALL 


= LE: 
nD ae ee ae eS a Se a Se ee a ene ees ee ee ees ee eo oe ee 











-~ t+ 
O/T; ieee Um 0.0 el.0 S1.5 1107 30.399 15.7 Se.4 


TIO /Ti’ 50.0 1.55 0.60 0.70 0.92 0.43 1.43 0.70 1.07 R 


TiO /Ti' 0.01 7.05 29.6 5.40 10.7 1.97 2.16 4.00 8.63 F 








TiO,/Ti* G7) =). = = 0.04 2. Zone 0.04 R 
BiOs/Ti me - = - 0.36 6 0.26 F 








Tio'/o. 286. 7.80 1.20 3.20 1.80 3.78 4.64 4.47 2.78 R 


TiO /O_ 0.07 35.3 58.5 26.3 20BeP 17-2 weuGn0G" 925.0 e 22-59 


ae GR Ca CS Ge ee GD en ee ee ee Gee ee 2 eee oe 





Ti0,/O” Oo. - = = 0.08 Bon) 0.09 R 
Ti0,/O° 1.60 = = = 0.82 = 0.66 ~ 0.71 F 


= =e, ee eG Sr Ge ee ae ee ne Oe ee 


NOTE: All quantities are in percentages. Experimental data 
taken from Reference 22. "R" and "F" refer to 
whether the mechanism involved was recombination or 
fragmentation. The TFA+B entries were obtained by 
averaging the data from the TFA and TFB positions. 
The ALL entires represent the average from all three 
adatom sites. 
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TABLE 23 


Comparison of Experimental and Simulation Results on 
Multimer Formation Mechanisms for Vanadium and Niobium 

















a Qe Can G-e-e-e eee ee eee ee ei ee a a ee a ee a ie 9 ne 6 | ee on eee SS ee 


NOTE: 


All quantities are in percentages. 


se &: 
YU (2x2) (2x2) 
eee STORE ATP 
10.0 10.8 Za 
PAT ST | 2.44 ae 
Gz02 7276 22a0 
OA = - 
0.78 - - 
DS. 220 iia 
0.19 90.4 76.0 
G24 - = 
Vee) - - 


P C 
RATIO YU (2x2) (2x2) 
mone SLOP ATOP __ 
-,,+ 
0” /Nb UNK 13.6 30.8 
++ 
NbO"/Nb” - UNK 4.55 5.01 R 
NbO /Nb* ——*UNK 6.04 14.3 F 
NbO,/Nb* UNK = 0.42 R 
NbO;/Nb UNK = 0.28 F 
NbO"/O UNK 32.7 16.6 R 
Nbo /O 0.72 43.3 47.4 F 
NbO,/0 UNK aa Te 8 
Nbo,,/O 14.0 = 0.88 F 


taken from Reference 23. 
the mechanism involved was recombination or 


fragmentation. 
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Experimental data 
"F" refer to whether 








Simulation. The experimental data was obtained by Yu 
[Refs. 22,23] using 500 ev neon and argon to bombard the 
three metals. His experiments used cesiated-oxygenated 
polycrystals and through the change in the work function 
induced by the cesium, he was able to obtain information 
on the charge state dependence of the two multimer ejection 
mechanisms. The ion ratios presented in the tables were 
calculated by extrapolating the relative yield of each 
Species back to a zero change in the work function which 
should represent the cesium-free case (Appendix J). 

No exact comparison can be made due to the uncertainty 
with respect to relative ionization probability, oxygen 
placement and density, as well as differences in ion mass 
and energy. However, if the ratios of ions formed by the 
two mechanisms were similar, the computer simulations would 
Support Yu's proposed model. This was not found to be the 
case. In each of the oxygen positions and densities the 
ratio of fragmentation to recombination was found to be 
Opposite that reported by Yu. 

This variance between the simulation and experimental 
results is probably not due to differences in lon energy 
since Table 12 indicates that the total multimer yield 
decreases as the ion energy increases. The percentage of 
multimers formed by fragmentation also decreases, but this 
is believed due to the longer-ranged transport of energy 


which allows more ejection mechanisms to operate. The 
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activation of more mechanisms could then decrease the rela- 
tive fraction of multimers formed by fragmentation. The 
simulation yield obtained at 500 eV was, in any case, too 
small for meaningful comparison. 

Failure of the simulations to support Yu's theory may 
be caused by an inadequate treatment of the ionization proba- 
bility contribution to the determination of multimer sta- 
bility. Since the program does not determine charge state, 
some of the multimers that were classified as stable may, 
in reality, be unstable and thus reduce the percentage 


found to be formed by fragmentation. 
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V. CONCLUSIONS 


The program has the inherent flexiblity to permit de- 
tailed analysis of virtually any aspect of Sputtering. The 
ease with which parameters may be isolated and varied 
greatly facilitates the determination of the significant 
factors involved, as well as the interdependence of these 
imaetOrs . 

The results of these simulations support earlier find- 
ings concerning the effects of channeling and that the majority 
of sputtered atoms originate in the top layers of the crystal. 
The angular distribution of ejected atoms agrees with ex- 
perimental results in the case of vanadium and niobium, 
however the discrepancies noted in the ejection patterns for 
Eitanzum requare further study. 

The decrease in oxygen yield from crystals with adatoms 
in the three-fold positions indicates that the oxygens are 
most stable in these sites. Such positions would also 
facilitate their incorporation into the metal lattice as the 
oxygen exposure is increased, resulting in the formation of 
a stable oxide in the top layers. It is also possible that 
the oxygens assume an A-top position at low oxygen exposures 
and as the coverage is increased, oxygen-oxygen repulsion 
causes a migration to the three-fold sites. 

The disagreement between the simulation results and 


Yu's theory and findings on multimer formation need not be 
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due to a failure of his basic model, but rather to an incon- 
plete understanding and treatment of the chemistry involved 
after the atoms are ejected from the surface by the simulation 
model. On the other hand it is somewhat difficult to under- 
stand how the yield of TiO’, presumably formed by recombination 
as suggested by Yu, can be so much larger than TiO , from 
fragmentation, considering the probabilities involved in two 
or more atoms finding themselves in the same location with 
relative velocities favorable enough to permit bonding. It 
may also be that a significant percentage of the negative 

ions stabilize themselves by electron ejection prior to 

their collection in the mass spectrometer. Yu has also 
pointed out that several discrepancies in his findings may 

be due to the fact that as the oxygen density is increased 
Tio’ and ites may be formed by both mechanisms. 

An item that bears further study is the increase in the 
most probable energy of the ejected substrate atoms when 
oxygen is introduced in the A-top position. At present 
there is no apparent reason for such an increase. If this 
were a model artifact one would expect the effect to carry 
over, to some degree, to crystals with adatoms in the bridge 
positions. Instead, the average and most probable energies 
for the bridge and three-fold sites are essentially the 


same. 
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APPENDIX A 


POTENTIAL FUNCTION GRAPHS AND PARAMETERS 
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Fig. 15. Ti-Ti solid phase potential function and associated parameters. 
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V-V solid phase potential function and associated parameters. 
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Fig. 17. WNb-Nb solid phase potential function and associated parameters. 
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Fig. 18. Ti-O solid phase potential functionand associated parameters. 
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Fig. 19. Ti-Ti gas phase potential function and associated parameters. 
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Ar-Ti potential function and associated parameters. 
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